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Abstract: In this work a lattice formulation of a super symmetric theory is proposed 
and tested that preserves the complete supersymmetry on the lattice. The results of a 
onedimensional nonperturbative simulation show the realization of the full supersymmetry 
and the correct continuum limit of the theory. It is proven here that the violation of 
supersymmetry due to the absence of the Leibniz rule on the lattice can be amended only 
with a nonlocal derivative and nonlocal interaction term. The fermion doubling problem 
is also discussed, which leads to another important source of supersymmetry breaking on 
the lattice. This problem is also solved with a nonlocal realization. 
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1. Introduction 

Supersymmetric theories have drawn much attention over the last decades. There is a 
variety of application for this symmetry ranging from supersymmetric quantum mechan- 
ical systems to extensions of the standard model (like MSSM) and string theory. While 
perturbative calculations and some nonperturbative arguments already relveal some in- 
vesting features of supersymmetric theories, it is still desirable to have a nonperturbative 
tool at hand for further investigations. Lattice calculations have proven to provide an im- 
portant insight in the nonperturbative sector of quantum field theories like QCD. To get 
reliable results of supersymmetric models with this method the theory must be discretised 
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in such a way that the correct theory is reproduced in the continuum limit. The estab- 
lished techniques of discretisation lead, however, to a breaking of supersymmetry on the 
lattice. There are various attempts to get, nevertheless, a supersymmetric theory in the 
continuum limit (cf. Q [2], Q for reviews). These approaches contain a fine-tuning towards 
the correct continuum limit [f|] and a partial realization of supersymmetry on the lattice 
H 0- However, the violation of supersymmetry at a finite lattice spacing seems to be 
unavoidable. The reasons for this violation of the continuum supersymmetry on the lattice 
are further investigated in this work. As a possible solution nonlocal realizations will be 
investigated. 

The most difficult problem is the absence of the Leibniz rule that is necessary for 
the invariance of the (interacting) lattice action. This problem was first discussed in ||, 
where as a possible solution an action with a nonlocal interaction term was presented. 
This approach was further pursued in (9|, 10 1, where also the breaking of the translational 
invariance of this suggestion was found. In [11] a discretisation of the four-dimensional 
Wess-Zumino model was suggested that respects all the continuum supersymmetry. Its 
disadvantage is a nonlocal derivative operator in addition to the nonlocal interaction term 
in the lattice action. In [^] a "No-Go theorem" was presented, that for the implementa- 
tion of full supersymmetry on the lattice either a nonlocal interaction term or a nonlocal 
derivative is necessary. According to this argument the nonlocal SLAC derivative should 
provide an realization with intact supersymmetry on the lattice. Such an action was al- 
ready investigated in [13], where despite the usage of the nonlocal derivative operator still 
a violation of supersymmetry at a finite lattice spacing was found. I will present here a 
more complete version of the simple "No-Go" statement. I will show that only with a 
nonlocal derivative and a nonlocal interaction term supersymmetry can be realized on the 
lattice. The conditions for a controlled approximation of the continuum theory on the 
lattice, presented here, also exclude the suggestion of Q. 

For the evaluation of the importance of the supersymmetry breaking by all local lat- 
tice realization, I will compare this effect with other sources of supersymmetry breaking 
on the lattice. These other important sources of supersymmetry breaking can be traced 
back to the fermion doubling problem. Although the considerations are quite generic for 
supersymmetric lattice realizations, I will focus here on Wess-Zumino type models. 

According to the "No-Go" statement the two options for the lattice simulations are 
that either locality or supersymmetry is not realized at a finite lattice spacing and must be 
recovered in the continuum limit. The best way would be to compare the results of both 
of these options. If they agree a supersymmetric and local continuum limit can be ensured 
in the nonperturbative sector theory. 

It demands, however, a large numerical effort to simulate the nonlocal supersymmetric 
lattice actions. I will present here a lattice formulation that allows for a simulation of 
at least the lowerdimensional models. The first results of a such a fully supersymmetric 
lattice simulation are shown in this paper. 

The paper is organized as follows: First I will give a short summary of supersymmetric 
models in the continuum, that are the starting point of the general discussion. In section ^ I 
will introduce a lattice discretisation for supersymmetric theories and show how a violation 
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of supersymmetry appears due to the absence of the Leibniz rule on the lattice. In section || 
a simple "No-Go" statement for local super symmetric lattice realizations will be presented. 
A concrete lattice realization that allows for a simulation with intact supersymmetry on 
the lattice will be considered in section ||. In most of the earlier work the supersymmetry 
breaking due to either the Leibniz rule or the fermion doubling problem was discussed. 
For a more complete discussion and a comparison of these effects, section [5] adds some 
considerations of the doubling problem. Section || contains some perturbative arguments 
for the local continuum limit of the nonlocal lattice actions. The more important nonper- 
turbative results of a supersymmetric lattice simulation are presented in section |7j The 
conclusions ^ at the end of the paper contain some suggestions for further improvements 
of the full supersymmetric lattice simulations. This work is based on an updated version 



of the results in [14]. 



1.1 Continuum supersymmetry 

The results of this work are quite generic for supersymmetric field theories, i. e. field theories 
with a symmetry that connects fermions and bosons and contain derivative operators in 
the symmetry transformations. Gauge theories demand some further discussions, especially 
when nonlocal operators are considered. Therefore, the main focus of the paper are Wess- 
Zumino type models. Models of that type are used for the matter sector of supersymmetric 
extensions of the standard model. The simplest example is the one-dimensional Wess- 
Zumino model. A derivation of the action in superspace and its off-shell representation is 
shown in appendix |A[ The on-shell action of this model is 

S =fdt Q(d t <^) 2 + $U + \(W'{^)) 2 + jjW"(<p)i/j\ . (1.1) 

It contains the scalar bosonic field (p and a complex fermion field tp. The bosonic poten- 
tial, the Yukawa interaction between fermions and bosons, and possible mass terms are 
derived from the superpotential W(tp) that is a polynomial of the fields. 1 Hence, there are 
relations between the bosonic and fermionic vertices. A possible mass parameter appears 
in certain bosonic interaction vertices. The action transforms under the supersymmetry 
transformations 

6<p = eil> + $e, 6il> = (dt<p-W'(<p))e, 5$ = -£(d t <p + W f (ip)), (1.2) 
into the integration of a total derivative term 

5S = -e[dt [(d t 1;)W'((p)+1>W , '(<p)d t <p] (1.3) 
''dt {(dS)W'(y) + ipW"(ip)d t ip] . (1.4) 



For the invariance of the action appropriate (e. g. periodic) boundary conditions must be 
assumed. In addition also the Leibniz (chain) rule is used. In the off-shell formulation of 
the theory the same terms appear as variation of the action. 2 

1 The prime denotes a derivative of the superpotential (W'(ip(t)) = ™W(p(t))). 

2 A contribution Ftp™ + mj)ip n ~ 1 ip transforms into (<9tV0v?" + n.'^ip n ~ 1 dtip since the auxiliary field F 
generates a derivative of a fermionic field and ip a derivative of a bosonic field. 
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There are similar models in other dimensions (with or without extended supersymme- 
try). The most prominent example for this class of models is the four-dimensional J\f = 1 
Wess-Zumino model. In the off-shell representation it reads ( fi~5|| ) 



S= (fx 



\d»Ad»A + \d„Bd»B + - \f 2 - \g 2 

+ m ^FA + GB- 
+ g (FA 2 - FB 2 + 2GAB - $(A - 



:i.5) 



This model contains the scalar field A, the pseudoscalar field B, a Majorana-spinor ip and 
the auxiliary fields F and G. 3 An important feature of super symmetric theories is the 
cancellation between bosonic and fermionic quantum correction. Bosonic loop diagrams 
are canceled by fermionic ones in perturbation theory. In the four-dimensional case this 
leads to the well-known nonrenormalisation theorem: To get finite expressions only a wave 
function renormalization (Z) that is the same for all fields is needed. The corresponding 
counterterm is 



^counter — J <fx(Z — 1) 



\d^Ad»A + \d^Bd»B + + X -F 2 + X -G 2 



(1.6) 



As in the one-dimensional case the action (L5) is invariant up to terms that can be identified 



with surface terms when the Leibniz rule is applied. 

More generally one can show that a nontrivial Lagrangian can be invariant under 
all supersymmetry transformations only up to a total derivative term. This is a direct 
consequence of the algebra: The commutator of some supersymmetry transformations 
applied to a field leads to a derivative of the field. All the transformation of a field can 
hence only be zero if it is constant. Schematically the action of a more general model can 
be written as 

^ =jd D x Qx(-d M cP + m 2 ) x + $(0 + m)il> 

+ mVi(x) + V 2 (x) + ^y{x)A , (1.7) 



The action is written here with the interaction vertices mV\(x) that carry in contrast to 
V2(x) besides dependence on the coupling constants also the mass parameter. The vertices 
of the Yukawa type interaction are summarized by Vy(x)- These vertices can involve a 
nontrivial dependence in the spinor indices. The x now can stand for several bosonic (e. g. 
scalar and pseudoscalar) fields. 4 



3 The problems of the correct definition of Majorana spinors in Euclidian space and on the lattice will 
not be discussed here. 

4 The condensed notation means x(~ d^d^ + m 2 )x = Y2r Xr(— d^d^ + m 2 )xr with, e. g., in the four 
dimensions M = 1 model the scalar field Yi = A and the pseudoscalar field X2 = B. 
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Supersymmetry transformations of such a general model lead to a variation of the 
action consisting of terms of the following form 



5S ote d D x \(d^) ■ ■ ■ (f) {ni) + ^(ft^ 2 )) • • • </> (nf) + ... + </> (1) • • • {d„4> {nt) ) \ ■ (1.8) 



0(1) . . . (J)( n i) denotes a product of fields, of which some may be of the same kind. One of 
these fields must be fermionic. As in the one-dimensional case each of these contributions 
vanishes if the Leibniz rule holds and appropriate boundary conditions are applied. The 
variation of the off-shell action are of a similar form. Just as in one dimension one can 
easily show in the off-shell formulation that (|l.8| ) comes from the variation of the terms 
with a product of rtf fields in the action. 

The appearance of these contributions can also be understood from a possible super- 
space representation of the theory. The supercharges contain Gramann and space-time 
derivatives (Q oc <% + Q^d^ with depending on the Gramann coordinates 9 a , and <% 
represents a derivative operator with respect to some of these coordinates). The relevant 
part of the action can be represented by a product of superfields #i, 02 3 • • •)) in- 

tegrated over superspace. Hence the supercharges generate the following variation of the 
action 5 

6S K / n d0 *d D x UeQ^)^ ■ ■ • $(™ f ) + $W(eQ*( 2 )) • • • ^ + . . .1 . (1.9) 

J a 

The contributions of the Grassmann derivatives in the supercharge vanishes after the Grass- 
mann integrations. The above expression (|1.8| ) can be identified with a component form of 
this variation. 6 

This brief review of continuum supersymmetry indicates that the generic form of the 
variation of the action is (|1.8| ). If the action should be invariant on the lattice as well as 
in the continuum we have to impose, as can always be done, periodic boundary conditions 
on the lattice. The second, more significant, condition for supersymmetry on the lattice is 
the Leibniz rule. 



2. The failure of the Leibniz rule and supersymmetry breaking on the 
lattice 

The crucial source of supersymmetry breaking on the lattice is the failure of the Leibniz 
rule for any discretised derivative operator. I will discuss possible solutions of this problem 
in this section. In the present discussion of the Leibniz rule the following assumptions are 
made: The same lattice derivative operator V^ 1 replaces the continuum derivative in the 
fermionic and bosonic kinetic part of the action as well as in the supersymmetry transfor- 
mations. As a starting point, the lattice action is chosen to be the same as the continuum 

5 Note that is written rather schematically in general one would have to take into account possible 
constraints that are applied to the superfields and several different supercharges. 

6 In other words: The variation of the highest component (highest power of Grassmann coordinates) of 
a product of superfields transforms into terms of the form presented in (|L.S|) as found in §1. 
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action ( |l,7|) except for a summation instead of the integration and the lattice derivative 
operator replacing the continuum derivative (the conventions of the lattice formulation 
are summarized in appendix |B|) . This derivative operator is assumed to be antisymmetric 
= —(V fJ- ) T and translational invariant. Its antisymmetric form ensures the hermiticity 
of the fermionic part of the action. The doubling problem of such an operator is discussed 
in section ||, since it introduces an additional source for supersymmetry breaking on the 
lattice. In that way, I want to separate the effect of the Leibniz rule from the fermion 
doubling problem. 

As detailed in the last section the Leibniz rule is needed for the invariance of the con- 
tinuum lattice action. When the continuum derivative is replaced by a discrete derivative 
operator the Leibniz rule becomes 

(V^ 1 ^))™ = </>«(V^ 2 )) m + (V^ 1 ))™^) , (2.1) 

for arbitrary lattice fields <^>W and <f>( 2 h One can easily verify the following statement. 



Lemma 1 The Leibniz rule, as presented in equation \2. \) , is violated by all lattice deriva- 
tive operators. 

Although this relation does not hold for any lattice derivative operator, it is still valid in 
the zero momentum sector, i. e. after the summation of m. 8 

Since the same derivative operator is used in the transformations the variation of the 
lattice action is just the lattice counterpart of ( |1.8[) , namely 



(nt)_ 



4 1} (V^ (2) )n • • • + ... + <t>$--- (V^>) J . (2.2) 



The fulfillment of the Leibniz rule after the summations ensures hence the invariance of 
a quadratic lattice action (rif = 2). In the higher than quadratic case the variation of 
the action (|2.2| ) is non-zero as for a product of three ore more fields the analog of ([^1]) is 
violated even after the summation. To show this in detail, let us consider such a product 



of fields in momentum space. The repeated application of (2J.) results in 



(V(^.--^))(p a ) = 

£ 5 L (p fcl + ...+p fcnf -p s )(V^(p fcl ) + ... + V^p fenf ))^ 1 )--.^). (2.3) 



7 To proof of this statement just insert a lattice field that is zero everywhere except on a single lattice 
point for </>"' and 4> ■ I thank A. Wipf for this simple argument. 



8 If one splits the antisymmetric matrix V according to (B.10) one can find the modified Leibniz rule 
for each of the components, J2 n VmUvi+rVn ') = (J2 n V™^)^^ + rf'fEn Vml^i 2 ')- Thus the 
Leibniz rule is fulfilled up to translations, which are not relevant under the summation. 

9 Cf. appendix M for the transformation to momentum space. 
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For the invariance of the action this equation must hold for p s = and all fields <fft' . This 
is equivalent to 10 

rif— 1 nf— 1 

^(Pfe + • • • +Pk nf w( Pkl ) + ... + v> fcnf )) = VM feJ - VM (E p^) = °- ( 2 - 4 ) 

i=l i=l 

The locality of the action is analyzed in the thermodynamic limit, where the discrete 
momentum becomes continuous, p ki —* Pi- 

For convenience the modulus of Yli=i Pi * s ^ rs ^ assume d to be smaller than the lattice 
cutoff A^ = for all directions fi. These assumptions fix the solution of equation fl2.4| ) 
to 11 

V*(p) = cu/* . (2.5) 

This derivative operator is (apart from irrelevant constants) the SLAC-derivative. The only 
possible solution is hence known to be nonlocal. This is the one of the basic observations 



of [Oj. It does, however, not mean that the Leibniz rule (2.1) is fulfilled for such a 
derivative operator. The solution is no longer valid, when the modulus of ^ST* Pi becomes 
larger than the lattice cutoff. The periodic continuation and the antisymmetry of in 
momentum space leads to a violation of ( |2.4D , when the argument of V^(p) is outside the 
first Brillouin zone. The violation is for (21^ - 1)A£ < Y^=i Pk z < + X ) A l ( V /^ lfl e Z ) 
given by 

rif— 1 rif— 1 

E^-E V ^) = 2A L ^^. (2.6) 

i=i i=i ti 

Even a nonlocal derivative operator is hence not enough to ensure the invariance of the 
lattice action. The basic reason of the violation of the Leibniz rule by the SLAC derivative 
is an incompatibility with the standard multiplication of the algebra of lattice fields. In 
this kind of multiplication a momentum of a product of fields, that exceeds the lattice 
cutoff is mapped back onto the region below this cutoff. It induces a periodic lattice "delta 
function" in momentum space (cf. ( |B.6D ), which is resolved by a periodic continuation of 
V(p). 

One can try to avoid this nonlocality of the lattice action with the introduction of a 
modified interaction term. For example, the product of three fields is represented on the 
lattice according to 

fc^)(x)^)(x)^)(x)° nth ^ tticc £ C mmmz ^M ■ (2-7) 

mi,m2,m3 

In Fourier space this ansatz brings equation ( |2.4D into the form (cf. ( |B.12| )) 

C(pi,p2,P3)(V fi {pi) + ^(p 2 )+V^(p 3 )) = Vm, (2.8) 



10 In momentum space V^ n is periodically continued for momenta larger than the lattice cutoff (cf. section 
11 ci = 1 in the continuum limit; additional constant contributions are zero because V M is antisymmetric. 
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which is solved by 

C(PUP2,P3) = II <5(V>i) + V> 2 ) + V*(P3)) • (2.9) 

This modification corresponds to a modified product rule on the lattice. This kind of prod- 
uct is defined in such a way that after the summation of the lattice index J} 5(^"^ 1 V M (pj)) 
instead of the usual (periodic) lattice delta Sl(Y17=~1 Pi) appears in momentum space. From 
the mathematical point of view the mentioned incompatibility of the Leibniz rule with the 
usual product on the lattice suggests such a modification of the product. 

For the symmetric derivative = (V^)^ such a solution was proposed in ||. How- 
ever, with a derivative different from ( |2.5| ) the solution breaks for continuous momenta the 
translational invariance on the lattice. A possible way out is to accept a modification of this 
invariance |J. Here this property is, however, considered to be even more important than 
locality. Another problem of the solution ( |2.9| ) with a symmetric derivative appears when 
discrete momenta of a finite lattice instead of the continuous momenta in the thermody- 
namic limit are considered. Then it corresponds to a projection onto the trivial solutions, 
Pjtj = 0, pk 2 = 0, or pk 3 = (as already mentioned in g). Consequently, the nonlocal 
interaction term allows for every finite lattice only a trivial interaction with at least one 
field at zero momentum, i. e. 

' ( p x ^{x)4> {2 \x)^{x) ° nt ™- £ [^m {2) {-p k )^\Pk) 

k 

+^ {1 H-^)0 {2) (O)0 {3) (^) + {1) (-p fc )^ {2) (^)^ (3) (O)] • (2.10) 

It is obvious that such an expression is supersymmetric. The fields with zero momentum 
arguments are similar to a mass term for the other fields in the product. Although for a 
continuous momentum the correct classical (tree level) continuum limit has been shown in 
Q, such a solution can not be useful for the lattice simulations. 12 The nonlocal interaction 
term alone does not provide a sufficient lattice realization of the theory. 

A translational invariant choice for a modified lattice product in ( j2.7| ) is C nitn2>n3 = 
CV ni _ n2 w m _ n3 ) (cf. ( |B.15| )). The resulting product of three fields then leads to 

£ On 3 -™- m2 [(vV 3) )m 3 0S0S + 0S(v^ (1) )m^S + ^ffi(vv (2) ) w 

-p 2 )C(pi,p 2 )(V(pi) + V(P2) " V(pi +P2))<P {1) (pi)^ 2) (P2) , 



mi,m2,m3 



Pl,P2,P3 

for the supersymmetry variation of the action. Obviously, this approach can be generalized 
to the situation of nf fields, and the condition for the invariance of the action, cf. ([2.4|), 
becomes 

"nf— 1 nf— 1 

£v(pi)-V(5» 

.1=1 1=1 



C(pi, . . . ,Pn f -l) 



0. (2.11) 



2 Perhaps this interaction term can, however, be interpreted as the interaction with a kind of mean field. 
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It corresponds to the Leibniz rule for rtf — 1 fields with the modified lattice product 

(0(D * ^ * ... * : = Y: C l ^..., l . mn{ _ 1 ^\ • • • tftgl 

mi,...,m„ f _i 

= £ C(,m 1 ,..,m„ rl ^ ) 1 • • • ^1 • (2.12) 

mi,..,ra„ f _i 

This time the modification of the product is defined such that after the summation of 
the lattice index the usual (periodic) lattice delta P<) i s combined with an 

additional factor C(pi, . . . ,p nf ) in momentum space. The modified product is commutative 
as C is assumed to be symmetric under the exchange of its arguments. 

With the help of the modified product on the lattice, we arrive at a more general 
realization of a supersymmetric lattice action. Although the Leibiz rule (|2.1| ) can not be 
fulfilled by any lattice derivative operator the generalized ansatz allows the construction 
of a supersymmetric realization of a nontrivial lattice action, as detailed later on. Before 
I present a concrete realization of such a supersymmetric lattice action, I first address 



the question in as much locality can be realized in such an approach. In [12] a "No-Go 
theorem" for a local supersymmetric lattice realization was found. Since the effect of the 
violation of the Leibniz rule by the SLAC derivative for large momenta was not included 
in these investigations, I get here an even stronger statement for the violation of locality 
by supersymmetric lattice actions. 

3. A No-Go theorem 

With the results of the discussion presented in the last section I now introduce the assump- 
tions for of a simple No-Go statement. 

As above the same antisymmetric lattice derivative operator is assumed to replace the 
continuum derivative in the action and the supersymmetry transformations. In addition 
to the lattice derivative operator a modified translational invariant product on the lattice 



( 2.12 ) is included. The modification of the product is assumed to appear in every term in 
the action, except for the kinetic part that contains the terms with the derivatives and the 
squares of the auxiliary fields. It appears also in the nonlinear part of the supersymmetry 
transformations. 

A kind of "smoothness condition" for the limit of infinitely many lattice points is 



assumed. This condition should exclude trivial solution similar to ( 2.10| ). When the number 



of lattice points increases, also the number of modes {^ nt \pk n{ )) that are allowed to 
interact with all modes of the other fields satisfying the constraint pk, l{ + YlT^ Pki = 
should increase. Otherwise it is hard to extrapolate the continuum limit, where no other 
constraint besides p n{ + Pi = i s imposed on the interacting lattice fields. In the 

thermodynamic limit at least in a certain interval around p n{ = all momenta should 
be allowed to interact as in the continuum without any further constraint. As the lattice 
spacing tends towards zero this interval has to increase. 

The locality of a lattice derivative operator (V^) nm means that it decays exponentially 
when the distance between the lattice points n and m increases. The width of this expo- 



-9- 



nential decay is proportional to the lattice spacing. A generalization leads to the locality 
condition for C . In this case locality means that the elements of Ci—rni,...,i—m n . 

decrease 

exponentially when the distance between the lattice point I and one of the mi, . . . , m n{ 
increases (at fixed distance between I and all the other mj). 

The careful definition of the conditions for a suitable lattice realization is the crucial 
part of the derivation of the following No-Go statement. Its proof is rather trivial. 



Lemma 2 (No-Go theorem) In order to get an interacting supersymmetric lattice the- 
ory one needs a nonlocal derivative operator and a nonlocal interaction term. 



The basic message of this statement is that either locality or supersymmetry must be given 
up on the lattice. 

For the proof of this statement I first recall some well-known facts that relate properties 
of a function with its Fourier transformed. If a function is differentiable to all orders, 
its Fourier transform decays faster than any polynomial. For an exponential decay it 
must be analytic in a region around the real axis. For the Fourier transformed of the 
lattice operators (in thermodynamic limit) these conditions must not only hold for the 
whole Brillouin zone. In addition they must be fulfilled at the boundary of the Brillouin 
zone, when the function is periodically continued. Hence locality means nothing but the 
analyticity of the periodically continued functions V^(p) and C(p\, . . . ,p nf -i)- In the case 
of C(pi, . . . ,p nf -i) this must hold with respect to each of the pi, . . . ,Pn f -l holding all 
other arguments fixed at an arbitrary value. If an analytic function is constant in a certain 
region (and not only on isolated points) it is constant everywhere. The same is true for 
the derivative of an analytic function. This implies that if the derivative operator has 
the momentum space representation (2.5) in a certain region of the BZ, it has to obey this 
momentum space relation in the complete BZ. Otherwise it would violate analyticity inside 
the BZ. It is hence a SLAC type derivative, that is analytic inside the BZ but nonlocal 
because periodicity is not fulfilled. 

The invariance of the action under supersymmetry transformations demands that equa- 
tion ( |2.11| ) is fulfilled in the interacting case, i. e. rif > 2. Therefore either the term in square 
brackets or the function C(p\, . . . ,p n{ -\) has to be zero. According to the smoothness con- 
dition there must always be a continuous region for the sum of the nf—1 momenta where 
C(p\, . . . ,p nf -i) is nonzero. In this region the term in the square bracket must be zero and 
the derivative has to be (2J5). A derivative operator that has this kind of momentum space 
representation in a certain region of the BZ can, however, not be analytic and periodic. 
Thus the derivative operator must be nonlocal. 

On the other hand, the term in square brackets can only be zero in the region where 
the sum of the nf — 1 momenta stays inside the first BZ. In this case the derivative operator 
is the nonlocal SLAC derivative. Hence at least in the region where the sum gets larger 
than the lattice cutoff C(p\, . . . ,p n{ -{) must be zero. Since it is zero in a certain region 
of its arguments it can not be analytic. Thus also the modified lattice product must be 
nonlocal. 
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In contrast to [12| the smoothness condition and the effect of the periodic lattic delta 
are considered here. This implies that neither a local derivative operator nor a nonlocal 
product alone can provide a sufficient supersymmetric lattice realization. 

A possible approaches to circumvent this No-Go statement is a fine tuning of the 
bare parameters towards the supersymmetric continuum limit. Instead of the complete 
continuum supersymmetry one can also realize only a part of the symmetry on the lattice. 
It might be that this reduces the necessary fine tuning, but it can lead to other problems 
(cf. [pi ). Like the Ginsparg- Wilson relation for chiral symmetry, a similar relation could 
lead to a solution for the lattice realization of the symmetry. Such a symmetry relation 



has been found for supersymmetry in [17|. However, a solution of this relation in terms of 



an interacting local lattice action was so far not constructed. In 12] a certain realization 
using matrices instead of fields was proposed. Although it allows to circumvent the No-Go 
statement it induces an infinite number of degrees of freedom in the continuum limit. 

Different from these attempts I will construct here a nonlocal lattice representation. 
This allows for a way to find a supersymmetric lattice realization. Instead of supersym- 
metry one can investigate how locality is restored in the continuum limit. At least in low 
dimensions no fine tuning is necessary for this restoration. I will show in perturbative 
calculations and numerical simulations that correct results can be obtained from such a 
nonlocal lattice theory. 

4. A realization of the complete supersymmetric on the lattice 

In this section I will construct a lattice action that preserves all supersymmetries. A 
similar lattice realization was constructed in [|ll[] in the thermodynamic limit starting from 
a momentum space representation of the theory. Here I will formulate the lattice theory 
in a different way, that allows for a computation in the simulations. 

I start with the off-shell representation of the theory. Since anyway a nonlocal lattice 
derivative is needed one can, for convenience, choose the SLAC derivative (V M = (\7 sljAc ^^ 
in the lattice action and the transformations. The translation of a product of continuum 
fields into a nonlocal lattice term, cf. ( 2.12| ), is done for each term of the action separately 



according to the number rif of fields appearing in the product. A quadratic term needs 
no further modification. For a number of rif > 2 fields (no matter if they are bosonic, 
fermionic, or auxiliary) the appropriate C(p kl ,p k2 , ■ ■ ■ >Pkn -i) mus t be included. The Cs 
can be decomposed into a product of their one-dimensional counterparts C\, 

c(p kl ,Pk 2 ,- ■ ■ ,Pkn f -x) = n^K^i • • • >pL f -i) • ( 4<1 ) 

As found above in a realization with the SLAC derivative the nonlocal product has to 
exclude all momenta with a sum that exceeds the lattice cutoff. The following C solves 
(2.11) and leads to a supersymmetric lattice action: 



^l{P kl ,P k2 ,---,P k ) — \ . • <4.^ 

1 1 otherwise 
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In such a way a nonperiodic delta is introduced in the lattice theory, 

£ (nf-1) (Pifci + p k2 + ■■■ + Pk nf ) ■= S L (Pk! + ■■■ +Pk nf )C(Pk 1 ,Pk 2 , ■ ■ • ,ft„ r i) • (4.3) 

This nonperiodic delta function is nonzero only if the sum of the nf — 1 momenta is equal 
to Pk n ■ A nonperiodic delta was also introduced in [fUjl . 

For the purpose of the lattice simulation, an appropriate representation of the modified 
product with a nonperiodic delta is needed. Since C can in higher dimension be represented 
as a product of its one-dimensional counterparts, it is enough to find a solution the one- 
dimensional case. Consider a lattice with the (nf— l)-fold number of lattice points compared 
to the original lattice. The lattice spacing of this finer lattice is a /{rif — 1). Hence the 
periodicity of a Fourier space delta of this lattice (<5 nf _1 ) is (nf— 1)A^. The boundary of the 
first BZ of this finer lattice cannot be reached with the sum of the momenta p kl ■ ■ ■ Pfe nf _! • 
They are hence never folded back and the 5 ni ~ 1 (p kl + ... + Pk, l{ ) is equal to one only 
when Yl^^i 1 Ph exactly matches p nf , and otherwise zero. Thus it represents a nonperiodic 
delta. A one-dimensional interaction term with a nonlocal product of nf fields can hence 
be represented on the lattice in the following way (cf. equation ( B.20| )) 

dx^\x) • • • on th ^ attice £ 5^-\ Pkl + ... +p kni )^(p kl ) ■ ■ ■ jWfaJ 

(rif-l)iV-l 



f 



ki,...,k nf m=0 



- £ e-^ m ^ + --- +p ^U il Hp kl )---^H Pk J. (4.4) 



It is clear that a representation of a nonperiodic delta is also provided when 5 m_1 with 
m > ri{ instead of <5 nf_1 is used. Consequently for all products <5 m can be used, where m 
the highest appearing power of the fields (nf). 

Going back from momentum to real space representation one arrives at 13 



on the lattice 



(n f -l)7V-l 

£ 4® ■ ■ • with = Y, 4® ■ (4-5) 



rif — 1 

n=0 m 

The ((rtf — l)iV) x N matrix J-^} translates all of the fields into fields on the finer lattice. 
It comprises a Fourier transformation on the lattice with N lattice points and an inverse 
Fourier transformation on the larger lattice with (rif — 1)N lattice points. The matrix 
elements read explicitly 

~j n .\ sin(7r(m — nlint — 1)) . ,, , H . . , „. 

= AT '-^-i >J -rr for m i=- n /(n f - 1 and 1 otherwise. 4.6) 

In this representation the nonlocal product is defined by 

(rif-l)JV-l 

C — ■ ■ • F( nf ) (A 7) 



rif — 1 

1 n=0 



d The ordinary definition of the Fourier representation of the fields </>l^ is employed 
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This approach can easily be generalized to the higher-dimensional case. One obtains 
a fully super symmetric lattice action, which can be used in numerical simulations. The 
result of the simulations can be found in the section ^. 

For the on-shell representation of the theory the approach can also be applied. The 
modification of the product must then be chosen according to the power of fields that 
appear in the corresponding variation of the action ( |1.8[ ). The same modification of the 
product must also be applied for the nonlinear part of the supersymmetry transformations. 

The considered lattice action can also be represented completely on the finer lattice 
with rif — 1 lattice points, if rif is the largest number of fields that appear in a product. It 
leads to a theory where the momentum of all fields is constrained below the cutoff A/(nf — 1) 
instead of A. 

According to the No-Go theorem there are also other possibilities to realize an in- 
tact supersymmetry on the lattice. In each of these representations the nonlocal product 
constraints the momentum in the interaction terms to a region, where the Leibniz rule is 
fulfilled. For a high enough momenta of the fields there are effectively no interaction terms 
present in the theory. The interacting theory hence effectively lives on a lattice with a 
larger lattice spacing. The advantage of the approach presented in this section is that a 
contribution of the noninteracting modes is not present. 

5. The standard discretisation and supersymmetry breaking on the lattice 
due to the fermion doubling problem 

In the last sections I have discussed the violation of supersymmetry due to the absence of 
the Leibniz rule. In this discussions I have assumed an antisymmetric derivative operator 
in the lattice action. It is well-known that such a derivative operator introduces a doubling 
problem. I will discuss the doubling problem here, separately. This will also explain the 
assumption of the same derivative operator for fermions and bosons. Let us therefore 
allow different derivative operators and masses for fermions and bosons. The continuum 
derivative operators in the fermionic and d^d^ in the bosonic sector are replaced by 
lattice difference operators Vnm and O nm . The continuum action fll.7| ) is represented on 
the lattice by 

Son =J<Px Q<K-n + m l)<t> + j>(ft + 

+m fe Vi(0) + V 2 (0) + i>V Y (0)V) • (5.1) 

In the simplest discretisation in the bosonic sector is the forward and backward derivative 
(□ = Ylfi^ii- i > ( m b)mn = m^mn)- The discretisation in the fermionic sector is more 
difficult. It is well-known that a naive discretisation, i. e. with V^, introduces a doubling 
of the fermion species in the continuum limit. In accordance with the Nielson-Ninomiya 
theorem |18| , [l9| l all local representations with an antisymmetric derivative operator share 
this problem. The only way out is an additional contribution that is symmetric in mo- 
mentum space (m^ w \—p) = m^ w \p)). The common way to add such a contribution to 
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the action is with the identity in spinor space. Then it resembles a momentum dependent 
mass. To remove the doubling problem, this mass diverges at the additional zero modes 
(doublers) in the continuum limit. The doublers have then no dynamic contribution to the 
action and are "freezed out" . The most prominent example of such a term is the Wilson 
mass. It is 



In the case of supersymmetry the doubling modes are problematic. In the naive lattice 
formulation they appear only in the fermionic sector. Consequently, there are no longer the 
same number of degrees of freedom in the bosonic and fermionic sector, and supersymme- 
try is broken. As described above one can try to resolve this problem with a Wilson-type 
mass term in the fermionic sector and set mj = m + rrS w \ Due to the additional mass 
term mS w ^ and the different derivative operators in the bosonic and fermionic sector su- 
persymmetry is, however, violated even in the free theory. Classically it is restored in 
the continuum limit, where the difference between the derivative operators and the addi- 
tional mass term vanishes. In the quantum theory divergencies of the loop contributions 
alter this classical behavior and violate supersymmetry even in the continuum limit. Thus 
neither the difference between different derivative operators for fermions and bosons nor 
the additional Wilson mass term can in principle be assumed to have no influence on the 
continuum limit. 

Indeed, the cancellations of bosonic and fermionic loop contributions are based on the 
relations between the bosonic and fermionic vertices and propagators. These relations are, 
however, not valid in this discretisations, and the continuum cancellation do not appear in 
lattice perturbation theory. They can be recovered only in the continuum limit. To find 
out the relevance of certain contributions in lattice perturbation theory for the continuum 
limit, one can use the lattice degree of divergence defined by Reisz [Q. 14 A one loop 
diagram in lattice perturbation theory is represented as an integral of the loop momentum p 
constrained to the first BZ. When the lattice degree of divergence is negative, the continuum 
limit can be obtained in a naive way: The integrand is taken in the limit of a — > and 
afterwards an unconstrained momentum integration yields the correct continuum limit. In 
such a naive continuum limit the contribution of the Wilson mass and of the difference 
between the derivative operators disappears (the essential part for this limit comes from 
the vicinity of p = 0). A typical one loop diagram with fermionic and bosonic contributions 



As usual for lattice perturbation theory the calculation is carried out in the thermodynamic limit. 
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that cancel in the corresponding continuum expressions is 




m f {p) 

~-E m v?V)v£ s V) + ("v(p)) 2 , 

m b (U - VjL s) V{f } + mj - ml) 
(-□ + m 6 2 )(-£ M vMv? ) +mJ) 

mb — nif 



(5.3) 
(5.4) 



The last term (5.4) is nonzero due to the Wilson mass term for the fermions. The bosonic 
mass parameter mj, that appear there comes from one of the vertices mV\. The lattice de- 
gree of divergence of this term is D — 1. Already in one dimension the necessary cancellation 
between fermionic and bosonic loops is hence not recovered in the continuum limit. This ef- 



fect was investigated in detail in plT (cf. also [22, [Lq]). It explains the non-supersymmetric 



continuum limit of a lattice realization with a straight forward implementation and a Wil- 
son mass only for the fermions. 

The first term (^^) is due to the difference between the bosonic and fermionic propa- 
gators. Its lattice degree of divergence, D — 2, is smaller. In one dimension it is not relevant, 
whereas in two dimensions it can lead to a supersymmetry breaking in the continuum. 

To avoid this kind of supersymmetry breaking it is necessary to use the same propa- 
gators for fermions and bosons (□ = (Vf)^(Vf) At ; mj = nib)- The derivative operator then 
introduces a doubling problem also in the bosonic sector that is removed with the same 
Wilson mass as for the fermions. It is crucial that this mass parameter appears also in the 
bosonic vertices of the theory (m^Vi = m/Vi). This nontrivial modification of the theory 
leads to the correct cancellations of the considered fermionic and bosonic loop diagrams. 
The Wilson mass is hence a modification of the usual mass term in the superpotential. 



Such kind of lattice formulations were already investigated, e. g., in [22, g^j. However, 
the additional vertices introduced by the Wilson mass can modify or violate symmetries 
of the model on the lattice as found in fl6|| . The SLAC derivative opens the possibility to 
circumvent this problem without a modification of the superpotential. Since it is nonlocal, 
it does not introduce any doublers, and a Wilson mass term is not needed. 

The described adjustment of the fermionic and bosonic mass terms and derivatives 
alone does not imply a supersymmetric continuum limit. This can be seen even from the 
perturbative point of view. The cancellation of fermionic and bosonic loop diagrams is 
present, as it is in the continuum. One can show that the nonrenormalization theorems 
of the superpotential in the off-shell theory still holds on the lattice as well as in the 
continuum. Nevertheless in the four-dimensional model the wave function renormalization 
is not the same for all fields [pj|| , even in the continuum limit. In this way the violation 
of the Leibniz rule is present in the lattice perturbation theory. Although the masses and 
derivatives are the same for fermions and bosons non-supersymmetric counterterms are 
needed to get a supersymmetric continuum limit. 
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6. Lattice perturbation theory of the nonlocal lattice theory 



In the last sections I have shown that only lattice realizations with nonlocal operators 
can be invariant under super symmetry on the lattice. An explicit supersymmetric lattice 
action can be constructed with such operators. In addition nonlocal operators do not need 
additional Wilson mass terms. Although these mass terms can be consistently included in 
the superpotential they still introduce a serious modification of the on-shell theory. 

On the other hand the locality of the continuum theory should be respected in the 
continuum limit. For the SLAC derivative in lattice QED Karsten and Smit have shown in 
1 24] that nonlocal and noncovariant counterterms are necessary to achieve this. 15 According 
to the usual argument this shows that nonlocal lattice actions are not allowed in the 
simulations. 

Let us first consider the SLAC derivative, but no modification of the product in the 
interaction terms. 16 As has been shown in p3| , the twodimensional Wess-Zumino models 
needs no nonlocal or noncovariant counterterms to achieve the correct continuum limit. The 
propagators of a theory with the SLAC derivative are the same as in the continuum. The 
only difference to the continuum perturbation theory is a lattice momentum conservation 
at each vertex (with a periodic delta) and the restriction of each momentum integration to 
the first BZ. As a lattice degree of divergence of the diagrams it is enough to take the degree 
of divergence of the corresponding continuum diagram. 17 The nonperiodic delta leads to 



so called "umklapp" contributions in the loop diagrams [25, ^j, where the momentum 



in one of these delta functions becomes larger than the lattice cutoff. In two dimension 
these contributions vanish in the continuum limit [Ol. For the situation in more than 



two dimensions it was claimed in [25| that that the correct continuum limit is obtained 
without any nonlocal or non-covariant counterterms and the "umklapp" contributions are 
irrelevant. 

The only difference of this model with the SLAC derivative and the full supersymmetric 
model of section |] is the absence of the "umklapp" contributions. All momenta that appear 
in the propagators are restriction to the first BZ due to the nonlocal product ( |2,12| ). Hence, 
apart from the restrictions of the integration of the loop momenta, the diagrams are the 
same as in the continuum. If the "umklapp" contributions are irrelevant the continuum 
limit is obtained with the same counterterms as in the case without a nonlocal interaction. 
This is certainly true for the one- and two-dimensional models. Recently a more explicit 
proof of this fact was presented in In addition it was pointed out that the constraints 
on the momentum integration lead to new additional counterterms in the four-dimensional 
theory. It is still not clear in as much these results are consistent with [25| and whether no 
nonlocal counterterms are necessary in the four-dimensional theory. 



J In 



a different conclusion was drawn that was criticized in [E61 



in equation (5.1). 



16 This means □ = V? , (m b ) nm = (m/) mn = mS mrl , Vf = V ! 
17 The degree of divergence defined in can, strictly speaking , not be applied here since it assumes 
that the integrands are periodic smooth functions of the momentum. 
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7. Simulation of a complete realization of supersymmetry on the lattice 
in a one-dimensional model 



The aim of the present work is not only to suggest a possible nonlocal realizations, as done 
in section |4|. The solution were also checked on the nonperturbative level. Therefore lattice 
simulations of the one-dimensional full supersymmetric model were performed. In this way 
I could measure directly signs of a completely realized supersymmetry on the lattice. The 



continuum model is chosen to be the same as in [13]. This will allow a direct comparison 
with the discretisations considered there. 

7.1 The discretised actions 

The starting point is a discretised version of the continuum action ( |1.1| ) 

S L = \ £ (( VSLA( V)n + W 'MD + £ & ( V - C + WMnm) ^ • (7.1) 

n n,m 

In the continuum theory the superpotential is chosen to be W(<f>) = -y^ 2 f V? 4 - The 
derivative operator is in all of the discretisations considered here the SLAC derivative. 
The results of a discretised version with different implementations of the Wilson mass term 
can be found in [O] . 

On the lattice the following supersymmetry transformations are considered 

5W<p n = e^ n ; J<V, = 0; 6^ n = -e(V SLAC <p + W' L {<p)) n (7.2) 
6®<p n = $ n e; 6^ n = (V SLAC V - W' L {tp)) n e ; S^^ n = 0. (7.3) 

The discretisation of the continuum transformations ( |1.2j ) is in this way chosen according 
to the assumptions of section |2| and || with the same derivative operator as in the action. 

In the simplest discretisation no modification of the lattice product is included. Then 
the interaction is represented by 

W' L {<p) n = m + g<p*, (7.4) 

and 

W'l(ip) nm = {m + 3gtp 2 m )5 nm . (7.5) 

This discretisation is called here the unimproved model. The action is not invariant under 
the two supersymmetry transformations. 

As in [13] a second discretisation of the model is constructed with the so called "Nico- 



lai improvement". The difference between ( [7.11 ) and the corresponding lattice action are 
discretisations of surface terms, which are assumed to vanish in the continuum limit. It 
reads 

5 " = \ E ( VSLAC V + W 'M)l + E ^ ( V - ° + W 'Mnm) V> m , (7.6) 
n n,m 

with the interaction vertices (|7.4j ) and ( |7.5[) as in the unimproved model. It is invariant 
under the supersymmetry transformation ( |7.2[) , but not under the second supersymmetry 
transformation ( |7.3[ ). This discretised lattice model is called here the improved model. 
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The masses of the improved SLAC model and the 
full supersymmetric model (m = 10, g = 100) 

19 | , , , , , r- 

18.5 - imp. SLAC fermions 
imp. SLAC bosons 
18 - full susy. bosons i — — > 

full susy. fermions i — r — ■ 

17.5 - 




15.5 - 
15 - 
14.5 - 



0.005 0.01 0.015 0.02 0.025 0.03 0.035 

lattice spacing a 

Figure 1: The masses for fermions and bosons of the full supersymmetric model in comparison 
to the results of the improved model presented in ]l3| . The predicted continuum value is 16.865. 
In addition the linear fit of the data of the improved model (fermions: 16.81 — 18.53a; bosons 
16.84 - 17.97a) is shown. 



Finally I present here the first simulations with a model that respects all of the super- 
symmetry on the lattice. According to the construction in section |I] the continuum bosonic 
interaction vertices are represented on the lattice with 

3AT-1 

w'M n = m(Pn + ^-Yl r£Uvm? , (7.7) 

m=0 

and the fermionic ones with 

3N-1 

W'l(<p) nrn = m5 nm + ag ^ ^mln^mlmi'fmi) 2 • (7.8) 

mi=0 

For the current choice of the superpotential the fields <p are 

= FniWm , (7.9) 
m 

with according to ((O]). With this choice for W' L the difference between the improved 
and unimproved discretisation vanishes already at finite lattice spacing. The model is called 
here the full supersymmetric model. 

7.2 The algorithm 

For the simulations the same algorithm as used in |ll|] was applied. Supersymmetry 
demands for the simulation of dynamical fermions. This leads to a contribution of the 
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fermionic determinant in addition to the bosonic part Sb of the action in the path integral 
measure, 

SeffM = S B (<p) - logdeb K f (<p) . (7.10) 



The HMC algorithm |2{J is applied. The update algorithm follows a molecular dynamics 
trajectory determined by 

on . on 

<£n = q , TTn = - q , (7.11) 



d7T n d(/7 n 



and the Hamiltonian 



N-l 



n=0 



The numerical solutions of the differential equations ( |7.11| ) are computed with a standard 
leap frog algorithm. The fermionic contribution dS e s(tp)/dtp n is calculated from 



J-(logdet^) = = * , (7.13) 

where 



<9(£ n dip n 

When W'l is diagonal, the trace collapses. This reduces the cost of the numerical calcu- 
lations with a diagonal interaction term. This is not the case for the realization with a 
nonlocal interaction term, where the trace includes a sum over all lattice points. Never- 
theless the simulation is still possible at least in lower dimensions. 

7.3 Numerical results: Masses and Ward-identities 

To get an indication of supersymmetry on the lattice the masses and Ward-identities of the 
theory are measured on the lattice in the same way as in |H| . The masses are obtained from 
the exponential decay of the fermionic and bosonic correlators. The fermionic correlator 
of the SLAC derivative shows appart from the exponential decay an additional oscillating 
contribution that vanishes in the continuum limit. To enhance the signal for the fermionic 
mass already at finite lattice spacing a filtering technique was applied (for details cf . |H| ) . 
In addition the following Ward-Identities were measured 



^n-m = {Vn^Hm) + ((<5 (1) <p n )$ m ) 

Ward identity 1: K^_ m = (^ m ) - (MV SLAC V) m ) - {<PnW L ((p) m ) (7.15) 
Ward identity 2: R%± m = - (<p n (V SLAC <^) m ) + {<p n W' L (<p) m ) . (7.16) 



If supersymmetry is realized both of them must be identical zero. 

The obtained masses show not much difference for the considered discretisations. The 
masses of the improved and unimproved model were already presented in [13]. A good 
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Ward-Identity 2 
Ward-Identity 1 



6 8 10 

lattice point n 



12 



14 



Figure 2: The results of a measurement of the Ward identity in the free theory at small lattice 
spacing. This identifies the deviation from zero due to the statistical or systematic errors. No larger 
deviation from zero than 0.0003 is observed. 



agreement between fermionic and bosonic masses was measured in both models. In com- 
parison to a discretisation with the Wilson derivative 18 the discretisations with the SLAC 
derivative showed a nearly perfect behavior. Already at finite lattice spacing the results 
were very close to the correct continuum value. 

Since the simulation of the full supersymmetric model demands a higher numerical 
effort, a sightly lower statistic was used than for the improved and unimproved model in 



1 13]. Nevertheless 10 5 to 4 x 10 5 independent configurations were obtained. The masses 
coincide with the improved model within the statistical errors as shown in figure |l|. As 
already observed in the earlier calculations the bosonic masses show a larger statistical 
error than the fermionic ones. Even the fermionic masses, that are obtained with a high 
precision, agree with the linear extrapolation of the masses of the improved SLAC model 
towards the continuum limit. Thus the correct continuum value is obtained with the 



nonlocal full supersymmetric lattice realization. In [13] it was found that realizations with 
the SLAC derivative show an almost perfect behaviour. Only a small deviation form the 
continuum value is observed at finite lattice spacing. This is also true for the model with 
the nonlocal interaction. 

As the masses show not much difference between the considered models and do not pro- 
vide a clear sign of supersymmetry breaking at a finite lattice spacing, the Ward-identities 



18 This means the symmetric derivative and the Wilson mass term in the fermionic and bosonic sector as 
discussed in section || 
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The Ward-idenities of the unimproved model 
(m = 10, g = 800, iV = 15) 
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Figure 3: The Ward identity of the unimproved model at a larger coupling. The deviation from 
zero is clearly above the errorbounds. 



are considered. At lower couplings and a smaller lattice spacing the Ward identities vanish 
in all cases. Therefore a smaller lattice and a large coupling was chosen. The configurations 
were obtained in eight independent runs with each 5 x 10 5 independent configurations. To 
estimate possible systematic errors the Ward identities of the free theories were measured, 
which are analytically zero. The result is shown in figure ||, and no larger deviation from 
zero than 0.0003 is observed. This implies that if a Ward identity stays within these bounds 
(illustrated by a shaded region in the plot), it can be assumed to be zero within the er- 
rorbounds. This is not the case in the unimproved model. If the coupling constant g is 
large enough a clear deviation from zero can be observed (figure |3|). In the improved model 
one of the two Ward identities is zero within the errors (figure ||). On the other hand the 
second Ward identity gets an enhanced deviation from zero. 

This indicates one of the problems of a partial realized supersymmetry on the lattice: 
In order to avoid the breaking of one of the supersymmetries a larger breaking of the 
remaining supersymmetries is generically induced. Note furthermore that the measured 
Ward identities of the unimproved model are related via Tiff = — For the fermionic 
part of the Ward identity this is due to ($ x ij)y) = — {ipyipx}- The bosonic part of the 
relation between the two Ward identities comes form the invariance of the path integral 
measure under the transformation t — > — t in the action. This symmetry is present in the 
continuum and the unimproved models. Thus a partial realization of supersymmetry can 
only be achieved at the cost of a breaking of another symmetry of the continuum model. 
Additional problems that can be induced by a partial realization of supersymmetry are 
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The Ward-idenities of the improved model 
(m = 10, g = 800, iV = 15) 
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Figure 4: The Ward identity of the improved model. The Ward-identity 1 is zero within the errors. 
On the other hand Ward identity 2 shows a deviation from zeronearly two orders of magnitude above 
the errorbounds. 



discussed in |flq] . 

A comparison of the different results obtained in [13] shows that for unimproved and im- 
proved models with a symmetric derivative and a Wilson mass in the bosonic and fermionic 
sector the same sign of supersymmetry breaking appears. 

In the full supersymmetric model both Ward identities are zero within the errors (figure 
These results verify that only in the full supersymmetric realization with a nonlocal 
interaction term supersymmetry can be realized at a finite lattice spacing. It provides the 
first precise measurement of a fully supersymmetric model on the lattice. 



8. Conclusions and outlook 

In this work I have shown that a realization of a complete supersymmetry on the lattice 
can only be achieved with nonlocal interactions and a nonlocal derivative. This fact is 
due to the breaking of the Leibniz rule by all lattice derivative operators. The nonlocal 
interaction term comes form a modified product on the lattice that restricts the momen- 
tum modes. A concrete implementation for the realization of the nonlocal product was 
presented. Another important breaking mechanism due to the fermion doubling problem 
is also not present in this nonlocal discretisation. A local continuum limit without nonlocal 
or noncovariant counterterms can be verified in lattice perturbation theory at least in lower 
dimensions. Despite the nonlocal form of the action simulations in supersymmetric quan- 
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The Ward-idenities of the full supersymmetric model 
(m = 10, g = 800, N = 15) 
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Figure 5: The Ward identity of the full supersymmetric model. Within the determined small 
errorbounds it is zero. 



turn mechanics were performed in this work. The realization of a complete supersymmetry 
on the lattice can be verified with a high precision. In this way the present work has proven 
that simulations with intact supersymmetry on the lattice are indeed possible. 

The general discussion of lattice perturbation theory in this work reveals a hierarchy of 
supersymmetry breaking mechanisms. Each of the breaking mechanisms becomes relevant 
when supersymmetry breaking counterterms are needed for a supersymmetric continuum 
limit. The first problem is due to the Wilson mass, that must be included in the fermionic 
sector for all local lattice derivative operators. The mass parameter is, however, not present 
in the bosonic sector. This problem is relevant already in one dimension. In the two 
dimensional case the differences between the fermionic and bosonic derivative operators 
(propagators) can also become a relevant source of supersymmetry breaking. Both of these 
effects can be avoided when the same derivative and mass operators are chosen for fermions 
and bosons. In the on-shell theory this requires also an unconventional modification of 
some bosonic vertices. In lattice perturbation theory the missing Leibniz rule becomes a 
relevant effect in the four-dimensional models. This makes a full supersymmetric realization 
attractive especially in this case, where, however, it still has to be proven that no nonlocal 
counterterms are needed in the continuum limit. 

On the nonperturbative level it is still not guaranteed that the violation of supersymme- 
try can be neglected even in lower dimensions. Thus a model with complete supersymmetry 
on the lattice is a good cross check for other simulations. If the simulations of a nonlocal 
model with intact lattice supersymmetry and a local model without lattice supersymmetry 
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leads to the same results, the recovery of locality and supersymmetry in the continuum 
limit can be ensured. 

On the technical side the simulations of the nonlocal lattice action require a large 
numerical effort. The formulation used in this work provides a possibility only in the 
lower dimensional case. Further improvements are, however, possible. The nonperiodic 
delta function ( |4.3| ) is realized via a periodic delta function on a finer lattice. One can 
also formulate the whole lattice theory on the finer lattice, where no modification of the 
product is necessary. On the finer lattice the only important requirement for the full 
supersymmetric theory is a restriction of all momentum modes below the cutoff of the 
coarser lattice. In addition one has to keep in mind that the theory effectively lives on 
a coarser lattice when the observables are computed. This opens another possibility for 
the calculation of the supersymmetric theory. The lower cutoff of the theory may be 
implemented via a momentum-truncation of the force, such that no higher momenta are 
generated in the update step. This must be ensured for (pseudo-) fermionic and bosonic 
fields. Instead of a sharp cutoff one can also try to realize a smooth version of it. Instead 
of the full supersymmetric model one can also try to find improved actions with a reduced 
violation of the Leibniz rule and a moderate suppression of higher momentum modes. 19 

A disadvantage of the present approach is that it can not be directly applied for 
supersymmetric gauge theories. The SLAC derivative requires a gauge invariant connection 
between distant lattice points and according to [24] nonlocal and noncovariant counterterms 
are needed in the continuum limit. Perhaps a gauge fixed version of this theory can be 
used in the simulations. 
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A. Superspace formulation and off-shell representation of the one dimen- 
sional theory 

The one-dimensional (bosonic) superfield has the following expansion in the Grassmann 
coordinates 9 9 

<Z>(t,9,9) = <p + 9~Tp + %b9 + 99F. (A.l) 

It contains the real bosonic field ip, the fermions ip and ip as well as the auxiliary field F. 
The supersymmetry transformations of the superfields are generated by the supercharges 

Q = idg + 9d t and Q = id e + 9d t . (A.2) 

19 A similar treatment was suggested in |30[. 
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This implies the following supersymmetry transformation for the component fields: 



Sip = ietp — iipe , Sip = (dtp — iF)e , 

8$ = e{d t p + iF) , SF = -ed t ip - d t $e . (A.3) 

A Lagrangian is constructed using the covariant derivatives 

D = id s -Odt, D = id e - Sdt , (A.4) 

that anticommute with Q and Q. The action is hence obtained from an integration over 
the whole superspace, 



dOdOdt 



h(t, 9, 9)K*(t, 9, 9) + iW(*(t, 9, 9)) 



(A.5) 

= jdt fyd t ipf - ijdtip + l -F 2 + iFW'(p) - i$W"(<p)ip\ , (A.6) 
where W($) (here W($) = m§ 2 /2 + 5$ 4 /4) is a polynomial in $ and K = \{DD - DD). 
B. Some conventions 

In the n direction there are assumed to be lattice points separated by the spacing a^. 
The number of lattice points is odd in each direction. If the index of N and a is not 
specified the same number of points and the same lattice spacing applies for all directions. 
Each lattice point is labelled by a vector n with D integer components running from to 

— The lattice point is defined as x n = Ylfi( n )^ a ^- is the value of the field at this 
point. The vector e M has zero components except a 1 in its fj, direction (x n +e M is the next 
neighboring lattic point of x n in jjl direction). 

The size of the lattice in jjl direction is, consequently, = N^a^ and its volume 
0,1, = Yl^L^. Periodic boundary conditions are assumed for this volume. If not further 
specified the sum Y^ n °f a lattice index is 

E-nk E (B.D 

n M \ (n)„=0 / 

The lattice index is labeled by n and m. n-i or rij are treated as individual lattice 
indices; only the (n) M or {n) v stands for a component of an index. Thus (ni) M is the 
component of n± in (i direction. The same applies for the dimensionless wave vector k. 

The derivative operators used here are defined as 

symmetric derivative (V^fV)n = 2^(^n+e M - V?n-e M ) 

forward derivative (V^ + V)n = ^-( ( ^n+e M - <Pn) , 

backward derivative (V^ ' p) n = —(tp n - p n - e ^ 

SLAC derivative (Vf^p) n = -V J^i-l)™"- 1 — - ^ 



sin(7r((n) M -0/A r M ) 
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The SLAC derivative 20 , [31, 32, 33], is derived from a discretisation of the Fourier space 
representation of the continuum derivative operator. 

Functions on the lattice can be represented in Fourier space according to, 

0n = J>G%K PfcX " (B.3) 
k 

^(Pfc) =J2^~ ipkXn , (B.4) 

n 

with (pk)fi = and the components of k are integers running form —(N^ — l)/2 to 

(Nfj, — l)/2. (pgX here stands for a scalar product of the two vectors.) If not further 
specified the above sum over k represents 

(fc) M = (7V-l)/2 
fc Af (&)„=-(iV-l)/2 

4>(pk) is periodic in p^, cp{pk) = <Xp*; + e^/27r/a^) V7 € Z and all directions The mo- 
mentum (pfc)u of the modes is inside the Brillouin zone (BZ) defined by BZ = {(pu)| |pJ < 

The Fourier transformation implies the following representation of the delta on the 
lattice 

h{x n - x m ) := £ e ^K— ) = W~%W 

fc A* 

fcCPfa -p te ) := ^e-^i-^ = n^M^)^!),^), > (B-6) 

n fj, 

where <5( n ) M m M is one f° r ( n )^ = "V mod and zero otherwise. The two delta functions 
are periodic: 5 L (p k ) = 5 L (p k + e^Tr/a^); 5 L {x k ) = 5 L (x k + e^lN^) VI £Z and fj,. 

For the lattice perturbation theory the thermodynamic of the lattice expressions is 
performed. Then one gets the following Fourier representation with the now continuous 
momentum p: 



d D p 

Hp) = Y,^ e ~ ipXn - (b 



HvY VXn ■■= / ttpy pxn (B.7) 

BZ (2vr)^ 



(n 6 Z D now runs over an infinite number of lattice points.) One can easily change 
between the Fourier representation on a finite lattice and in the thermodynamic limit. The 
expressions in Fourierspace remain the same, only the discrete momentum p k has to be 
replaced by the continuous p. This continuous momentum is, because of the finite lattice 
spacing still restricted to the BZ. Consequently the delta in Fourier space is still periodic 
as on the lattice. 



3 This derivative is also called DWY derivative (Drell, Weinstein, and Yankielowicz). 
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In a similar way the Fourier representation of translational invariant operators with 
two indices is derived on the lattice: 

V£ m =J2^(Pk)e ipkiXn - Xm) . (B.9) 

k 

The obtained operator has the same periodicity as the Fourier space representation of a 
field. On the lattice the matrix entries of V should be real. Therefore, the imaginary part 
of V(p) is antisymmetric (9V(— p) = — 9V(p)) and the real part is symmetric (3?V(— p) = 
»V(p)). 

An antisymmetric derivative operator acting in [i direction can always be represented 

as 

N-\ 

VL=E ( v( v(r) )m«» with (V {r) )% m = 5 n+re ^ m -6 n - re ^ m , (B.10) 

r=l 

and some constants c r . 

For the operators defined in flB.2|) the Fourier representation is thus 

i 



V£ s) (p) = — sin(^a M ) 



(Vi t +) V(-))(p) = ^sin 2 (^/2) 



"i "Hp) = ^2 — sm 2 (p^ afl /2) 



^ AC (p)=i P ,. (B.ll) 

Apart from the SLAC derivative these are all functions of the type a d F(ap) where d is 
determined by the dimension of the operator. Obviously the behavior of F in the vicinity 
of p = is the important part in the continuum limit unless there are other points with 
F = 0. 

Now consider a more complicated operator C miv .. m . Translational invariance implies 
that a shift of all lattice points by x m is irrelevant. One way for the representation in 
Fourier space is 

C(p kl ,...,p kni ) = Yl Snw** e-^ x ^ + - + ^i x ^ . (B.12) 

mi,...,m„ f 

The translational invariance means for the Fourier space representation 

C( Pkl ,. ■ ■ ,p knt y iPkl+ '" +Pknt)Xm = CW • • ,Pk n{ ) • (B.13) 
This implies p^ = —p kl — ■ ■ ■ — P k 1 mod for all \i. A representation that explicitly 
implies the translational invariance is obtained by going from C with rif indices to a matrix 
C with rif — 1 indices according to 

Cm 1 ,...,m n{ = C(m„ f -mi),...,(m„ f -m„ r i) • (B-14) 
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A representation in Fourier space of this matrix is 

C mi ,...,m nf ^ = £ C( Pkl ,. . . , Pfcnf _ i)e ^ 1 ^ 1 + - + ^ f - 1 ^ f - 1 . (B.15) 

Thus the two representations are in Fourierspace related by 

C(p kl ,- ■ ■ ,Pk nt ) = 5 L [p kl + • • • +Pk nf )C(p kl ,. ■ ■ tPk^-t) ■ (B.16) 

With these operators the modification of the product on the lattice (|2.12 ) is defined. 

The commutativity of the product is ensured by the fact that C is invariant under 
the exchange of its arguments. A requirement for a modified product, not discussed in the 
main text, is its associativity, 

* 00) * cj ) W) l = ((^W * 0( 2 )) * (B.17) 
In Fourier space associativity on the lattice demands 

C(pk 1 ,Pk2,Pk 3 ,Pki) = ^2C{p k5 ,Pk 2 ,Pk 3 )C(p kl ,-pk 5 ,Pk 4 ) (B.18) 

or C{p kl ,Pk 2 ,Pk 3 ) = J }2h(Pk 4 ~Pk 2 -Pfci)C , (pjfc 4 ,pfc 3 )C(pfc 1 ,pfc 2 ). (B.19) 

Ay 

This condition is fulfilled for the proposed lattice action in section |4|. 
The matrix J r ^ n, t> used in this construction follows from 

/ 2itk ( a 



= £ -p ( • 



am — -n 



aN ^ V aN \ (nt - 1) 

■ 1it(N-\) f __n_\ 

* 2iv ^ m n^-iJ _ g i2n-(m-n/(n/-l)) 
~~ ^/V 1 - e i27r/AT(m-n/(n / -l)) 

sin(7r(m — n/{nt — 1)) 
aN sm(TT / N (m — n/(nj — 1)) 

It is clear that for this kind of matrix the following summation rule holds 

nfN-l 



(B.20) 



^ _ ^nm 1 ^~rim 2 ~ °~( x m 1 %m 2 ) ■ (B.21) 

:d 

E^ (1) ^ (2) )' = E^M 2) - CB-22) 



n=0 

This implies for the modified product of two fields 

,<'■>, ^> 

jp( n f) maps the fields tp® with a lattice size N on the fields ipn = ^ n Fmn <Pn with 
a lattice size (n/ — l)iV. The fields on the larger lattice have, however, a momentum 
constraint that is constraint below ^ instead the larger lattice cutoff n a ■ Thus 
generates a one to one map of p^ (pk) onto the modes of the larger lattice below the cutoff 
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